Anomalous dipion invariant mass distribution of the Y(45 ) decays into Y(15 )n + n and T(25 )n + n 
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To solve the discrepancy between the experimental data on the partial widths and lineshapes of the dipion 
emission of Y(45 ) and the theoretical predictions, we suggest that there is an additional contribution, which 
was not taken into account in previous calculations. Noticing that the mass of T(45) is above the production 
threshold of BB, the contribution of the sequential process T(4S) — > BB — > Y(nS) + 5 — > T{nS ) + n + n~ 
(n = 1,2) may be sizable, and its interference with that from the direct production would be important. The 
goal of this work is to investigate if a sum of the two contributions with a relative phase indeed reproduces the 
data. Our numerical results on the partial widths and the lineshapes dT(T(45) — > T(25, lS)n + 7r~)/d(m n +^) are 
satisfactorily consistent with the measurements, thus the role of this mechanism is confirmed. Moreover, with 
the parameters obtained by fitting the data of the Belle and Babar collaborations, we predict the distributions 
rfT(T(45 ) — » T(25, 15 )n + n~)/d cos 6, which have not been measured yet. 
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The dipion emissions of the bottomonia were observed a 
long while ago, and as an important channel to investigate 
the decay mechanism and inner structure of bottomonia, they 
have undergone careful studies from both theoretical and ex- 
perimental sides. The significance of the transitions stimu- 
lated a series of theoretical efforts J3-Sl ( se e Reft. 0-01 f° r a 
review). Before Y(4S) joined the game, the transitions occur- 
ring among only the bottomonia Y(3S) and Y(15,25) were 
below the BB threshold 0]] . With T(45 ) being experimentally 
observed, the situation becomes more intriguing because the 
BB channel is open. It is similar to the case of ^(3770), which 
is a charmonium with a mass just above the DD threshold. 
The CLEO-c Collaboration explored the dipion emission of 
Y(4S), however, only upper limits for Y(4S) — > T{\S)n + n~ 
and T(45) -» T(2S)n + n- were obtained $U&. With the data 
accumulation of the two B-factories, more precise measure- 
ment on the dipion emission of Y(4S) became possible. 

In recent years, a great progress on measurements of 
the T(45) dipion transitions has been made UU-HH. The 
Belle Collaboration reported their observation of Y(4S) — » 
T(lS)7r + 7r EMU whose partial decay width is F(Y(4S) -> 
Y(lS)7r + 7r) = (3.65 ±0.67(stat)±0.65(syst))keV. The BaBar 
Collaboration also presented their measurements on Y(45) 
decays into T(lS)7T + n~ and H{2S)n + n~ lfl3ll . and the partial 
widths are respectively T(Y(45) -> Y(lS>r + 7r) = (1.8 + 0.4) 
keV and T(Y(4S) -> r{2S)n + n~) = (2.7 + 0.8) keV HI . 
Besides measuring the decay rates of Y(45) — > Y(15)7r + 7r" 
and Y(45) -> f{2S)n + n ~, the Belle and BaBar collaborations 
also obtained the dipion invariant mass distribution, which 
is defined as dT(T(4S) — > Y(2S, 15 , )7r + ^)/c/(m ff + ff -) versus 



m n + n - and dT(Y(4S) — > Y(25, 15')7r + 7r _ )/(icos 9 versus cosB, 
where 9 is the angle between Y(45 ) and ii~ in the ji + tt~ rest 
frame. 



As shown in Fig. 3 of Ref. 11311 . a comparison between 
the experimental data and theoretical prediction yt] on the 
m n+n - distributions for Y(45) -> Y(15)?r + ^ _ and Y(45) 
Y(25 )n + ji~ decays was given. The measured distribu- 
tion for Y(45) — > T(lS)7T + n~ is close to that predicted in 
Ref. H, 

whereas the measured lineshape of the m^+^- dis- 
tribution for the Y(45 ) — > Y(25 fully deviates from the 
prediction of the Kuang-Yan model |30- Considering that the 
model is completely successful in describing heavy quarko- 
nium transitions Y(25) -> Y(15)?r + ^, T(35) -> Y(15)^ + tt- 
and if/(2S) — > J/i[m + 7T~, one should conclude that there must 
be some mechanism, which has not been taken into account in 
the old calculations. 



Moreover, the ratio fi = 



T(Y(4S)- 
T(4^T(T(15)^^-) 
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^gg- 1.16 was given 
by the BaBar Collaboration Ill3ll . Generally speaking, the 
Y(15')^ + ^ _ channel possesses a larger phase space than that 
of Y(25')^ + 7r _ , and it should result in a ratio % to be smaller 
than unity. As a good example, the rule is well satisfied for 
Y(35) as the measured % = 0.577 111 411 . The experimental 
value of the ratio R for Y(45 ) obviously declines from the rule, 
which might also imply existence of some uncounted factors. 
The very recent experimental result on the ratio H of Y(55 ) is 
given by the Belle Collaboration as R — 1.44 lfl5ll . which also 
contradicts to the simple argument. Thus, motivated by these 
phenomena, one should conjecture that anomalous production 
rates and lineshapes of Y(5S,4S) -» Y(lS)£r + 7r, Y(2S)7r + 7T 
are due to the available open BB channels lfm[l7ll . 

Y(45) is obviously different from the lower states. First, 
the mass of Y(45) is just above the BB threshold, so that 
the coupled channel effect becomes important. Another fea- 
ture is that Y(4S ) predominantly decays into BB pair because 
the mode is not suppressed by the OZI rule. Thus, when ex- 
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ploring hidden-bottom decays of Y(4S), we cannot ignore the 
coupled-channel effect, which may play a crucial role in fact. 

In general, there exist two mechanisms for the dipion emis- 
sion 

Y(4S) -» T(rtS)(p!>r + <>2>r<>3) with (»=1,2). 

One is a direct process, where two gluons are successively 
emitted and eventually hadronized into two pions, which is 
depicted in Fig. [T](a) and can be described by the QCD Multi- 
pole Expansion (QME) model yfl. Applications of the method 
have been discussed in all details, thus in this letter, we do not 
repeat the calculation for the direct transition, but alternatively 
introduce an effective Lagrangian to take care of this contri- 
bution. 

The transition amplitude of the direct transition Y(4S ) — > 
Y(nS )n + 7T~ can be written as 



Direct 

<f{n) /r / 2m 2 \ 

= -JT^S) ■ ey { „ s) {[q 2 - K { "\AM) z (l + —fj 

■ 



S-« 



(1) 



which was suggested by Novikov and Shifman in studying 
ip' — » J/ifm + 7t~ lfl8ll and the subscripts (S-wave) and (D- 
wave) denote the S-wave and D-wave contributions respec- 
tively. The mass difference between T(4S) and T(nS) is ex- 
pressed as AM. q 2 = (p 2 + pi) 2 = nh. denotes the invariant 
mass of n + n ~, while 9 is the angle between Y{4S) and if in 
the n^if rest frame. The pion decay constant and mass are 
taken as f„ = 130 MeV and m„ = 140 MeV, respectively. In 
Eq. ([T), k and T are left as free parameters when fitting the 
experimental data. 
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FIG. 1: The diagrams relevant to T(45) hidden-bottom decay with 
dipion emission. Here, Fig.[T](a) represents T(45) direct decay into 
Y(nS)n + n~ (n = 1,2), while Fig. [2(b) denotes the intermediate 
hadronic loop contribution to T(4S) — > Y(nS)n + n~ . 

There also exists another mechanism shown in Fig. Q] (b) 
contributing to the Y(4S) — > Y(nS)7T + n~ transitions, where 
the intermediate hadronic loop constructed by B and B mesons 
performs as a bridge to connect the initial state Y(4S ) and final 
state Y(nS)7T + 7T~ . This manifests the coupled channel effect, 
an important non-perturbative QCD effect. In other words, 
Y{4S) — > Y(nS)n + ii~ can be approximately expressed as a 
sequential decay process. Y(45') first transits into Y(nS) and 
scalar state S, then S couples with dipion and the picture is de- 
picted in the right panel of Fig. Q] For Y(4S) -> Y(lS)7r + ;r 
and Y(4S) -» T(2S)n + n- processes, S are cr(600),/ (980) 



and S = {<x(600) respectively, which are of proper quantum 
numbers and allowed by the phase space of the correspond- 
ing modes. Since the intermediate B and B mesons can be 
on-shell, as usually considered, the contribution from the ab- 
sorptive part of the loop is larger than the dispersive part, thus 
let us just keep the imaginary part of the loop. This definitely 
brings up errors, but in our strategy, we fit the data to obtain 
the model parameters, so that part of the uncertainty would be 
compensated. By the Cutosky cutting rule, the transition am- 
plitude corresponding to the absorptive part of the hadronic 
loop is written as 



At[T(4S) - 

M 



BB 



327r 2 M T( 4S) 
xM[B + (ki)B-(k 2 ) 



Y(nS)S — > Y{nS)n + n Jpsl 
J dnW[r(4S) -> B + (k 1 )B-(k 2 )] 

Y(nS)( Pl )S(q)] 



q 2 - m 2 s + ims Ts 



= [gosVv(pi ■ q) + gtoP^y]- 



J 1 V 

T(4S) T(nS) 

— m\ + ims F 



-g Sm P2 ■ Pi, 



(2) 



where the factor £ = 4 is due to charge conjugation B =± B 
and isospin transformation B° =± B + and B° =± B . Using the 
effective Lagrangian approach, we get 3\*[Y{4S) — > B + B ] = 



^T(4S) 
^T(«S)Bb4(«S) 



T(nS)S] = 



(-iki p -ip p )T [m B , p~] with the monopole form 

factor T[m 2 B ,p 2 ~\ = (A 2 - m 2 B )/(A 2 - p 2 ), where mg is the 
mass of the exchanged meson. The variable A is usually re- 
parameterized as A = nig + chAqcd with Aqcd = 220 MeV. 
This monopole form factor is introduced not only to reflect 
the structure of the effective vertex for the B + B — > Y(nS)S 
scattering process, but also concerns the off-shell effect of the 
the exchanged Z?-mesons at the t-channel. The coupling con- 
stants for the Y(45 )BB is estimated by fitting the partial de- 
cay width and gY(4S)BB — 

24, while for Y(nS )BB (n = 1,2) 
and SBB, the coupling constants are directly taken from Ref. 
lfl9tl . In addition, the coupling constants between the interme- 
diate scalar meson and the final 7r + 7r~ are g am = 16.2 GeV 1 
and gf nn = 2.40 GeV 1 , which are determined by fitting the 
corresponding partial widths fH l20ll . Finally we parameter- 
ized the whole hadronic loop contribution as a simple and ex- 
plicit Lorentz structure presented in the last expression of Eq. 
10, where g^} and gQ, corresponding to S-wave and D-wave 
contribution respectively, are extracted from the hadronic loop 
calculation and dependent on the dipion invariant mass . 

With the above preparation illustrated above, the total tran- 
sition amplitude of Y(4S ) — » Y(nS ~)n^n~ is expressed as 



AfTotal=AW + 2 e! ^ )A1 ! 



FSV 



(3) 



where the phase </>^ 1 is introduced by hand to account for a 
possible relative phase between the direct transition and the 
hadronic loop contributions. As a three-body decay, the dif- 
ferential decay width for Y(4S) — » Y(nS)n + 7T~ reads as, 

dT = 7 — 4 \M tn vJ 2 dml,^^dml +n - (4) 



(27r) 3 32m\ {4s) 



l T(nS)n L 
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with m\ (n?)n+ = (pi + p 2 f and m 2 ^^ = (p 2 + pi) 2 . The over- 
lined bar indicates an average over the polarizations of Y(4S ) 
in the initial state and sum over the polarizations of T(nS ) in 
the final state. The parameters about the concerned resonances 
are listed in Table. U 



TABLE I: The resonance parameters adopted in this work [l|,|2j 





Mass (GeV) 




Mass (GeV) Width (GeV) 


Y(15) 
Y(2S) 
Y(45) 


9.469 
10.024 
10.579 


tr(600) 
/o(980) 


0.526 0.302 
0.98 0.07 



In our present scenario, we keep four parameters T w ', /e (1) , 
$r , 0y o ( 98O ) f ree f° r Y(4S) — » Y(lS)7T + 7r~ and three free pa- 
rameters T {2 \ k {1 \ <p {2) for T(45) -> Y(2S)7r + 7T . So far, there 
are three measurements on the dipion invariant mass spec- 
trum of Y(45) -> Y(lS)7r + 7T given by Belle flH and 
BaBar fl3ll also measured the m n + n - distribution of Y(4S) — > 
Y(2S);r + ;r~. Our strategy is following. We have the Belle 
and BaBar data on the partial widths of Y(4S ) -> Y(2S )7r + 7T 
and Y(45) — > Y(lS);r + 7r~, and the measure dipion distribu- 
tions of dT I dm^+jj- , then using our model we make a best fit 
to the lineshapes when keeping the partial width to be consis- 
tent with the measured values. Namely, we will check whether 
the Belle ITlllT2Tl and BaBar data ITOl can be reproduced by 
our model. Our numerical computation is done with the help 
of the MINUIT package. 
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FIG. 2: The dipion invariant mass distributions for Y(4S) — » 
Y(1S)7T + 7T and Y(45) -> Y(2S)7r + 7T transitions. Here, the red solid 
curves are the best fit to the experimental data (blue points with error 
bars) QHIi in our model. Here, the data points with error bars in 
(a) and (b) correspond to the measured results given by Belle in Ref. 
fllll and a published work fl3l . respectively. 



For illustrating how to fit the line shape of Y(45) — > 
Y(15, 25 step by step, let us provide some details about 
the changes of theoretically predicted line shapes (see Fig. 
[3} of Y(4S) — > Y(lS)7r + 7r~ while considering the interfer- 
ence shown in Eq. ([3J. Since both <t(600) and /o(980) con- 
tribute to Fig. [TJ (b), |AiTotail 2 can be decomposed into six 
terms \M D \ 2 , |/VU 2 , 2Re(\M* D M tT ), |M/ | 2 , 2Re(M^M fo ) and 



2Re(M* T Mf ), where we abbreviate Afoirect 
A^and M^l ->M A . 
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FIG. 3: As an example, we illustrate the changes of the line shape 
(black solid lines) for fitting the BaBar data (blue dots with er- 
ror) of Y(45) — > Y(lS)7r + 7r~ flDl by adding the contributions from 
all the six modulus-square of individual amplitudes (|A4dI 2 > IA-U 2 , 
2Re(\M^MA IA-1/ol 2 , 2Re(M' D M fo ) and 2i?e(M;A* /o )) one by one. 
Here, the red dashed lines correspond to these modulus-square of the 
individual amplitudes. 



It is noted that for presenting the dipion invariant mass dis- 
tribution, the experimental colleagues use a different quantity 
"entries/mass", which exactly corresponds to the distribution 
c/F(Y(45) -> T(2S, \S)n*ir)ld{m^ r ), where the different 
lines correspond to the parameters obtained by fitting differ- 
ent data bases as indicated in the figures. Thus when we plot 
the dependence of the partial width over m^-, we deliber- 
ately normalize our numbers for a real comparison with the 
available data. 



TABLE II: The values of parameters for the best fit to the experimen- 
tal invariant mass distributions of Y(45) — > Y(lS)7r + 7r~. Here, the 
values listed in the second, third and fourth columns correspond to 
the best fits shown in diagrams (a), (b) and (c) of Fig. [2] respectively. 
The last row lists the partial decay width for Y(45) — > Y(lS)7r + 7r~ 
obtained by using the central values of those parameters. 



Belle 



-fur 



Belle (new) 



BaBar 



0.045 
0.249 
4$ 3.124 

Cso) 1-207 
T(keV) 2.255 



0.002 0.063: 
0.041 0.311: 
1.156 -2.64: 
0.489 0.973: 



0.003 0.040 + 0.001 
0.018 0.269 + 0.023 
0.609 -3.14 ±0.324 
0.723 0.918 ±0.183 



3.690 



1.799 



As shown in Fig. [21 we carry out the best fit (red solid 
curves) to the experimental data of Y(45) — » Y(lS);r + 7r~ and 
Y(45 ) — > Y(25 )ji + jt~ . The central values along with the er- 
rors of the fitted parameters are listed in Tables. Ull and Bill 



As shown in Fig. [21 the Belle data presented in Ref. ifTHl 
indicate that a double-peak structure appears in the range of 
m n +„- = 0.8 ~ 1.2 GeV, while later measured results given in 
Refs. lfl2[[l3ll do not show an obvious double-peak structure. 
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TABLE III: The parameters obtained from the best fit to the dipion 
invariant mass distribution of Y(4S) — > H{2S)n + n^ . The last row 
gives the partial decay width for Y(4S) — » Y(2S)7r + 7r~ obtained in 
our model. 
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0.838 ± 0.042 k {1> 0.664 ± 0.029 
-2.923 ±0.141 T(keV) 2.642 
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FIG. 4: The predicted dT/d cos 6 of Y(4S) -> Y(lS>r + 7r (left panel) 
and Y(4S) — » Y(2S)7r + 7r~ (right panel) transitions with the obtained 
values of the fitting parameters listed in Tables. 1 1 1 i 1 1 1 1 



This difference may be caused by different bin widths of the 
experimental settings of the Belle ifTll [l2tl and BaBar lfl3ll . 
Thus, in this work we separately fit their different experimen- 
tal data. Further experimental measurement will test whether 
there exist this double-peak structure in the dipion invariant 
mass spectrum of T(4S) — > Y(lS)7r + ;r~. 

In addition, we notice that the experimental data of 
T(4S) -> Y(lS)7!- + 7r and Y(45) -> Y(2S)7r + 7r are quite 
apart from each other. Our analysis indicates that consid- 
ering background and intermediate <x(600) and /o(980) con- 
tributions to T(4S) — > Y(lS)n + n~, both of the experimen- 
tal data ifTIUll can be well met. For Y(4S) -> Y(2S)7r + 7r 
process, there exist both background and intermediate cr(600) 
contribution, and it is different from the case of Y(4S) — > 
T(\S)n + 7i~ . Thus, the interference between background 
and intermediate cr(600) contributions reasonably results in 
a smooth line shape of the distribution of the dipion invariant 
mass spectrum of Y(4S) — > Y(2S )7r + 7r~ . 

It is noted that with the parameters obtained by fit- 
ting the lineshapes of the Belle and BaBar collaborations 
on DYIdm n , n - for Y(45) -> Y(l S )7r + 7T and Y(45) -> 
Y(2S)n + 7T~, we predict the angular distribution dF/dcos 9 as 
shown in Fig. 3. Such distributions have not been measure 
yet, so that the prediction will be tested by the further mea- 
surement of the Belle group or even a more accurate measure- 
ment at the proposed Super-B factory. 



As aforementioned there are several puzzles for the di- 
pion emission of T(4S). The anomalous ratio of 'R = 
T(45)^r( Y( 1 5 5F?r] ~ 1 - 16; and the lineshapes dF(T(4S) — > 
Y{2S,lS)7T + 7T~)/dm„+„- cannot be explained by the QCD 
multipole expansion theory, even though the theory works 
well for obtaining dipion emission rates of lower Y states. 
Noticing the mass of Y(4S) is just above the production 
threshold of BB channel, it is natural to conjecture that on- 
shell BB is an intermediate state and the sequential process 
Y(45) -» BB -> T(nS ) + S -> Y(«5) + tt+tT may play a cru- 
cial role. Namely this sequential process may have a sizable 
amplitude and interfere with the direct dipion emission from 
T(4S). The sum of the the two contributions (the two ampli- 
tude may have a relative phase, which is fixed by fitting data 
in this work) indeed results in a satisfactory explanation to the 
data. From Fig. (2), one can note that our model satisfacto- 
rily reproduce the lineshapes given by experimental measure- 
ments. Moreover, in Fig. (3), our prediction on the angular 
distribution dF/dcosd will be testified by the further mea- 
surements of the Belle-II group and/or maybe, the proposed 
Super-B factory. 

This result confirms our conjecture that the long-distance 
effect, which can be treated as a final state interaction in- 
deed plays a crucial role. Basically the direct dipion emission 
occurs via the OZI suppressed process, whereas, the open- 
bottom process Y(45) — > BB does not suffer from the sup- 
pression. Thus even though the final state interaction is re- 
alized via hadronic loops, its amplitude may be comparable 
with the OZI suppressed direct dipion emission, so that the 
interference is important. For lower Y(nS) (n < 3), BB can- 
not be on-shell, the absorptive part of the hadronic loop does 
not exist and the contribution may be negligible. 

As a conclusion, the puzzles of the dipion emission of 
T(4S) can be satisfactorily removed by taking into account 
the final state interactions. 
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